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CONTRAST DUAL FREQUENCY IMAGING 

The present invention relates to ultrasound imaging techniques, and in 
particular to ultrasound imaging techniques that use contrast enhancement 
5 agents to enhance discrimination between different structures under 
examination within a target medium or body, e.g. tissues and other structures 
within the human body. 

Over the past decade, contrast agent in a form of tiny gas bubbles has been 
10 introduced to improve image quality in ultrasound imaging techniques. The 
gas bubbles are infused into a region of interest to increase backscattered 
echoes from desired organs of interest. To extend the utility of ultrasound 
contrast imaging,, research has been actively focused in developing efficacious 
-ultrasound contrast agents and new contrast- specific imaging techniques. 

Second harmonic-based techniques enhance the detection of ultrasound 
contrast agents within many structures, such as cardiac chambers. Second 
harmonic techniques exploit the differences in response of gas microbubbles 
and tissue to ultrasound irradiation. Soft tissues are known to be linear 
20 reflectors of ultrasound energy, whereas contrast bubbles exhibit a non-linear 
or harmonic behaviour when interacting with ultrasound waves. 

By coincidence, it was found that imaging tissue at the second harmonic 
frequency, even without injection of contrast bubbles, improves image quality 
25 considerably. Tins mode, called native tissue harmonic imaging, has 
demonstrated advantages in various clinical applications. 

However, recent studies have shown that imaging contrast bubbles and tissue at 
second harmonic frequency only is also associated with some disadvantages, 
30 such as harmonic contamination from tissue. Tissue harmonic contamination 
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reduces the contrast detection capability and therefore the contrast to tissue 
ratio. Higher harmonic frequencies were then suggested as a way to increase 
the ability to detect or discriminate between contrast gas bubbles and 
surrounding scattering medium and have shown promising in vitro and in vivo 
5 results, as described in Bouakaz et aL, Ultrasound in Medicine and Biology, 
Volume 28, Issue 1, pages 1-137 (January 2002). 

It as an object of the present invention to provide an ultrasound imaging 
technique that is able to improve contrast bubble detection over existing 
10 -techniques. 

According to one aspect, the present invention provides an apparatus for 
ultrasound imaging comprising: 
a signal generator; 

15 a transmit transducer coupled to the signal generator for simultaneously 

irradiating a target with a relatively low frequency ultrasound conditioning 
signal and a relatively high frequency ultrasound detection signal; 

a receive transducer for receiving echo signals from the irradiated target; 

and 

20 a signal processor adapted to process the received echo signals to detect 

the presence of first structures within the target causing a first magnitude of 
detection signal echo arising from periods when the conditioning signal is in a 
first phase and a second, different, magnitude of detection signal echo arising 
from periods when the conditioning signal is in a second phase. 

25 

According to another aspect, the present invention provides an apparatus for 
ultrasound imaging comprising: 
a signal generator; 
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a transmit transducer coupled to the signal generator for simultaneously 
irradiating a target with a relatively low frequency ultrasound conditioning 
signal and a relatively high frequency ultrasound detection signal; 

a receive transducer for receiving echo signals from the target; and 
5 a signal processor adapted to process the received echo signals to 

differentiate first structures within the target whose physical properties change 
as a function of the conditioning signal from second structures whose 
corresponding physical properties remain substantially invariant in response to 
the conditioning signal. 

10 

According to another aspect, the present invention provides a method of 
ultrasound imaging comprising the steps of: 

simultaneously irradiating a target with a relatively low frequency 
ultrasound conditioning signal and a relatively high frequency ultrasound 
15 detection signal; 

receiving echo signals from the target; and 

processing the received echo signals to detect the presence of first 
structures within the target causing a first magnitude of detection signal echo 
arising from periods when the conditioning signal is in a first phase and a 
20 second, different, magnitude of detection signal echo arising from periods 
when the conditioning signal is in a second phase. 

According to another aspect, the present invention provides a method of 
ultrasound imaging comprising the steps of: 
25 simultaneously irradiating a target with a relatively low frequency 

ultrasound conditioning signal and a relatively high frequency ultrasound 
detection signal; 

receiving echo signals from the target; and 

processing the received echo signals to differentiate first structures 
30 within the target whose physical properties change as a function of the 
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conditioning signal from second structures whose corresponding physical 
properties remain substantially invariant in response to the conditioning signal. 

Embodiments of the present invention will now be described by way of 
5 example and' with reference to the accompanying drawings in which: 

Figure 1 shows • variation in a contrast agent bubble radius under 
insonification with a conditioning signal (the radius-time curve), for a bubble 
having nominal (resting) radius of 2 microns; 

Figure 2 shows the echo signal response of a contrast agent bubble to a 
10 detection signal as a function of bubble size, for radius 1.85 microns (figure 2a) 
and radiixs 2.34 microns (figure 2b); 

Figure 3 shows the echo signal response of a contrast agent bubble to an 
insonifyixig ultrasound beam having conditioning signal frequency of 0.5 MHz 
and detection signal frequency of 3.5 MHz; 
15 Figure 4 shows the echo signal response of a contrast agent bubble to an 

insonifyimg ultrasound beam having conditioning signal frequency of 0.5 MHz 
and detection signal frequency of 3.5 MHz, the echo signal being filtered at the 
second harmonic of the detection signal frequency; 

Figure 5 shows an optical image generated to determine the resting 
20 radius of a contrast agent bubble; 

Figure 6 shows the variation in a contrast agent bubble diameter under 
insonification with a conditioning signal (the diameter-time curve) for the 
bubble of figure 5 having nominal (resting) diameter of 4 microns, as measured 
with an ultrafast camera; 
25 Figure 7 shows filtered diameter-time curves for the bubble of figure 5, 

filtered at 0.5 MHz and 3.5 MHz; 

Figure 8 shows a schematic diagram of exemplary insonifying 
apparatixs; 

Figure 9 shows first and second successive transmit waveforms used in a 
30 conventional pulse inversion imaging technique; 
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Figure 1 0 shows first and second successive transmit waveforms used in 
a pulse inversion conditioning and detection signal technique according to the 
present invention; 

Figure 1 1 shows first and second echo signals from non-linear structures 
5 irradiated with the waveforms of figure 9; 

Figure 12 shows first and second echo signals from non-linear structures 
irradiated with the waveforms of figure 1 0; 

Figure 13 shows the sum of the echo signals of figure 11; 
Figure 14 shows the sum of the echo signals of figure 12; and 
10 Figure 15 shows a schematic diagram of an echo phase shift correlation 

process. 

The method of the present invention is based on blending two frequency 
components for simultaneous transmission to irradiate or excite a region of 

15 interest of a target body under examination (generally described herein as a 
'target')- A first, relatively low frequency component signal, referred to herein 
as the 'conditioning' or modulation signal, is combined with a second, 
relatively high frequency component signal, referred to herein as the 'detection' 
signal. The lower frequency (conditioning) signal is used to modulate the size 

20 of contrast agent bubbles by making them vibrate. The high frequency signal, 
which is blended with the low frequency signal, is used to interrogate the 
modulated bubbles, e.g. to enable imaging or detection by way of the local 
oscillations of the contrast agent bubbles that are induced by the low frequency 
signal. 

25 

An advantage of mixing the conditioning signal and the detection signal resides 
in the fact that the. conditioning signal alters the physical properties of certain 
contrast agent structures (e.g. gas bubbles). An important physical property of 
the gas bubbles that is modulated by the conditioning signal is size, and 
30 therefore also resonant frequency. The size of a gas bubble therefore varies 
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according to the compression and rarefaction caused by the acoustic pressure 
applied by the conditioning signal. The size of a gas bubble will reach a 
minimum during the positive pressure cycle of the conditioning signal and a 
maximum during the negative pressure cycle of the conditioning signal. In 
5 other words, the bubble size will oscillate, expanding in the negative 
(rarefaction) phase of the conditioning signal and shrinking in the positive 
(compression) phase of the conditioning signal. 

During these two phases of the conditioning signal, the detection signal is 
10 transmitted and used to detect or image the same bubbles. Hence, the detection 
signal will sense the same bubbles but at two different stages: small and large 
size depending on the phase of the low frequency signal. During the expansion 
phase (rarefaction of the lower frequency signal), larger bubbles are insonified 
•with the higher frequenc3' detection signal, and during the shrinkage phase 
15 (compression of the lower frequency signal), smaller bubbles are insonified 
with the higher frequency detection signal. Consequently the response of the 
gas bubbles to the detection signal will be different in the positive and negative 
cycles of the conditioning signal. 

20 "When non-oscillating scatterers (such as body tissue) are present (no gas 
bubbles present in the target medium), the response will be identical in both 
phases of the conditioning signal since these scatterers do not oscillate. This 
finding will increase the distinction in response between gas bubbles and tissue, 
-usually termed contrast agent to tissue ratio compared to other existing 

25 methods. 

With reference to figure 8, to deploy such a method, an ultrasonic array 
transducer or 'probe' 80 is used that is able to transmit both the higher and 
lower frequency components. Such a requirement can be fulfilled with, for 
30 example, the dual frequency transducer described in Bouakaz et al, Ultrasound 
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in Medicine and Biology, Volume 28, Issue 1, pages 1-137 (January 2002). 
The dual frequency transducer described therein contains two different types of 
active elements that are sensitive at two distinct frequency bands. Other types 
of transducer probe may be used 5 as will be apparent to the person skilled in the 
5 art. 

The transducer 80 is connected to a signal generator comprising first and 
second signal sources 81, 82 and a signal blender 83 for summing or otherwise 
combining the relatively low frequency conditioning signal from source 81 and 
10 the relatively high frequency detection signal from source 82. The blended or 
combined signal is supplied to the transducer 80 which irradiates a region of 
interest or target 84 in a body 85 under examination. 

The target 84 includes different types of structure that are capable of being 
15 resolved and separately imaged using the effects of the transmitted ultrasound 
signals. The target preferably includes at least: (i) a first type of structure that 
has at least one physical property which varies as a function of the conditioning 
signal; and (ii) a second type of structure for which the corresponding physical 
property does not vary significantly as a function of the conditioning signal. ■ 

20 

In preferred applications, the first type of structure is a contrast agent, such as 
fluid or fluid-filled 'bubbles 5 , which produce a non-linear response to the 
ultrasonic excitation. Preferably, the bubbles are gas bubbles as well known in 
the art of ultrasound imaging. 

25 

In preferred applications, the second type of structure is body tissue which 

produces a substantially linear response to the ultrasonic excitation. It will be 

i. 

understood that the response of the first and second types of structure may be 
dependent upon the frequency and / or intensity of ultrasonic excitation, and 
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therefore suitable such parameters for the ultrasonic excitation may be selected 
accordingly. 

The contrast agent may be introduced into the target body using known 
5 techniques that cause the contrast agent to selectively locate at selected 
locations in the target. 

A receive transducer 86 is used to receive echo signals £ E* from the target 84, 
which echo signals are passed to a signal processor 87. The signal processor 87 

10 is adapted to process the received echo signals E to differentiate between first 
structures within the target whose physical properties change as a function of 
the conditioning signal and the second structures whose corresponding physical 
properties remain substantially invariant in response to the conditioning signal. 
Various techniques that may be used by the signal processor 87 will be 

1 5 discussed hereinafter. 

A control unit 88 may be provided to coordinate and synchronise the activities 
of the signal generator components 81, 82, 83 and the signal processor 88. 
Those skilled in the art will recognise that the functions of the transmit and 
20 receive transducers 80, 86 may be combined in a single transducer or 
transducer array. 

Preferably, the system includes an image display unit 89 adapted for the display 
of the first and / or second structures and their relative locations within the 
25 target body, according to known principles of ultrasound imaging. 
Alternatively, or in addition, the image display unit could include any suitable 
form of output device, including a data storage receptacle. 

As used in the present specification, the expressions 'imaging 5 and 'detection 5 
30 in the context of the transmitted and received ultrasound signals is intended to 
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encompass the obtaining of any kind of location data relating to specific 
structures within the target 84. 

Referring now to figure 1, the radius-time curve 10 is shown for a vibrating 
5 contrast agent bubble that is being insonified with 0.5 MHz signal having an 
acoustic amplitude of 180 lcPa. This low frequency signal is used as a 
conditioning signal (produced "by signal source 81, figure 8). The bubble has a 
resting or nominal radius of 2 microns as indicated by the flat portion 1 1 of the 
curve. The bubble is encapsulated by a shell whose elastic properties are 
10 similar to those of known, available contrast agents. The bubble expands and 
compresses under the influence of the conditioning signal. It reaches a 
maximum radius of 2.34 microns (reference numeral 12) during the expansion 
phase and shrinks to a minimum radius of 1.85 micron (reference numeral 13) 
during the compression phase. 

15 

A higher frequency detection signal, generated by signal source 82 is then used 
to image the oscillations of this bubble at its maximal and minimal size. Figure 
2(a) shows the response 20 of* the small (1.85 micron) bubble, and figure 2(b) 
shows the response 25 of the large bubble, to a 3.5 MHz detection signal with 
20 an acoustic amplitude of 120 kPa. Since the smallest bubble has a natural 
resonance frequency closer to the detection signal frequency (3.5 MHz), its 
acoustic response demonstrates much higher amplitude 21 than that of the 
larger bubble (shown at 26). In addition, the smaller bubble clearly 
demonstrates more non-linearity in the response. 

25 

Figure 3 shows the response 30 of the bubble having nominal 2 micron radius 
to an excitation signal that contains dual frequencies 0.5 MHz and 3.5 MHz 
(simultaneous transmission). The dotted line 31 represents the radius-time 
curve (shown as a function of change in radius) of the bubble to the signal 
30 having only a single frequency 0.5 MHz. This shows the expansion phase 
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(positive-going portion 32) and the compression phase (negative-going portion 
33) of this bubble relative to the nominal radius (zero portion 34). The solid 
line 35 shows the acoustic response at 3.5 MHz, which indicates a strong 
asymmetry. 

5 

The bubble exhibits a stronger response during the compression phase 33, 
which corresponds to the smaller bubble radius (down to 1.85 microns): 
During the expansion phase 31, where the original bubble grows to up to 2.34 
microns, the bubble response is weaker. Figure 3 demonstrates that the bubble 

10 response is different depending on the phase of the conditioning 0.5 MHz 
signal. During the compression phase of the 0.5 MHz signal, the bubble 
shrinks and gets smaller. At this time the detection signal at 3.5 MHz will 
detect or image a small bubble. At the rarefaction phase of the 0.5 MHz signal, 
the bubble expands and the detection signal will detect a different bubble with 

1 5 much larger size . 

Thus, in a general sense, it will be noted that the detection signal echoes have a 
first average magnitude (relatively large) arising from periods (e.g. 33) when 
the conditioning signal is in a first (compression) phase and a second average 
20 magnitude (relatively smaller) arising from periods (e.g. 32) when the 
conditioning signal is in a second (expansion) phase. 

Figure 4 shows the same results but with these, the echo signals 40 are filtered 
around the second harmonic frequency (7.0 MHz). The dotted line 41 
25 represents the radius-time curve (change in radius) indicating the response of 
the bubble to the 0.5 MHz signal. The solid line 42 shows the second harmonic 
(2H) echo signal. When the bubble is smaller and closer to resonance size 
(negative-going portion 43), it generates more non-linearities and thus more of 
the second harmonic component compared to when it expands (positive-going 
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portion 44) and the non-linear behaviour of the bubble is damped. The dotted 
line 41 is the same as line 31 in figure 3. 

Ultra-fast optical observations of bubbles insonified simultaneously at two 
5 different frequencies were carried out using an ultra-fast camera. To 
interrogate the bubbles at two different frequencies, two different transducers 
were used. The transducers had centre frequencies of 0.5 MHz and 3.5 MHz. 
The transducers were mounted in a Plexiglas container and positioned such that 
their focal distances coincided at a depth of 75 mm. The optical observations 
10 were earned out at a frame rate of 13.45 million frames per second and 128 
successive frames were recorded. The acoustic pressures as measured 
separately were approximately 100 kPa at both the 0.5 MHz and 3.5 MHz 
frequencies. Commercially available contrast agent was used. 

15 Figure 5 shows the first frame 50, which was recorded at the arrival of the 
ultrasound signal. Bubbles with various sizes were present in the view and the 
following analysis was based on the bubble 51 which had a resting diameter of 
approximately 4 microns. 

20 Figure 6 shows the global diameter-time curve 60 of that bubble 51, which was 
oscillating under the effects of both 0.5 MHz conditioning signal and 3.5 MHz 
detection signal. The curve 60 shows clearly two frequencies present in the 
bubble response, a low frequency indicated by period 61 and a higher 
frequency indicated by period 62. Figure 7 shows the filtered diameter-time 

25 curves 71, 72 around 0.5 MHz (dotted line 71) and 3.5 MHz (solid line 72) 
frequencies. 

Up to 3 microseconds (period indicated at 73), the 0.5 MHz signal amplitude is 
low and the vibrations of the bubble 51 dialing this period of time are weak. In 
30 addition, the compression and expansion phases of the bubble are small 
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meaning that the bubble size does not change significantly between the 
compression and rarefaction phases 74, 75. Consequently the vibrations 72 of 
the bubble at 3.5 MHz are almost equal between the two phases 74, 75. 
However beyond 3 microseconds, the 0.5 MHz signal amplitude has increased 
5 and much stronger vibrations of the bubble now occur at this frequency. The 
bubble 51 radius decreases to about 3 microns and expands up to 6 microns. 
During these oscillations, the 3.5 MHz detection signal interrogates the same 
bubble 51 of variable size. This is seen in the response 72 at 3.5 MHz, which 
shows different responses depending whether a compression phase (period 
10 shown at 76) or an expansion phase (period shown at 77) of the 0.5 MHz 
conditioning signal is seen by the bubble. 

It will be understood that selection of appropriate conditioning signal 
frequencies and intensities and selection of appropriate detection signal 
15 frequencies and intensities will be dependent upon the physical properties of 
the body under examination and of the first and second types of strucixires 
therein. 

In a typical application of imaging human or animal tissue, in which 
20 conventional, commercially available contrast agents are used, the preferred 
conditioning signal frequency lies in the range 10 kHz to 5 MHz. More 
preferably, the conditioning signal frequency lies in the range 100 kHz to 800 
kHz. The detection signal frequency preferably lies in the range 1 MHz to 50 
MHz, and more preferably lies in the range 1 MHz to 5 MHz. 

25 

The signal blender 83 for combining the conditioning and detection signals 
may use simple signal summation, mixing, amplitude or power modulation 
techniques. 
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Signal processing of the echo signals in the signal processor 87 to detect the 
differences between the positive and negative phases of the conditioning signal 
can be performed using a number of techniques. These include pulse inversion, 
Doppler or correlation methods as will now be described. 

5 

Pulse inversion, Doppler or correlation methods can all be applied to: (i) the 
full echo response; (ii) the second harmonic of the echo; (iii) higher harmonics 
or sub-harmonics of the echo. Depending upon the exact nature of the body 
under examination, different degrees of discrimination between the first and 
10 second types of structure within the target may be achieved. 

With reference to figures 9 to 14, the use of a pulse inversion technique for 
signal processing will now be described. Figure 9 shows first and second 
pulses 90, 91 for insonification of the target 84 according to known techniques. 

15 The target 84 is excited using a first transmit signal pulse 90 followed by a 
second transmit signal pulse 91 that is inverted (180 degree phase shifted). In 
figure 9, the first and second pulses each comprise an envelope of effectively a 
single frequency. In the present invention, the same pulse inversion sequence 
is used, but in this case transmit pulses 100, 101 essentially each comprise a 

20 blended signal comprising two frequencies: the conditioning signal and the 
detection signal. 



Figure 1 1 shows first and second echo signal pulses 110, 111 resulting from the 
first and second excitation signal pulses 90, 91. Figure 12 shows the first and 

25 second echo signal pulses 120, 121 resulting from the first and second blended 
excitation signal pulses 100, 101. It is immediately clear that the difference 
between the first and second echo signals 120, 121 is substantially larger than 
the difference between the first and second echo signals 110, 111, illustrating 
the substantial improvement in discrimination possible with the present 

30 invention. 
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The echoes from the pulse inversion process are summed according to known 
techniques and the results shown in figures 13 and 14. 

Figure 13 shows the echo summation (difference signal) 130 for the 
conventional process, indicating only small differences between the first and 
second (inverted) echo signals. This difference represents second (and/or 
higher harmonic components that are not cancelled by the summation process. 
In case of a fully linear system, the summation would be zero. 

Figure 14 shows the analogous result 140 for the double frequency imaging 
(conditioning signal frequency plus detection signal frequency) technique of 
the present invention. The echo difference signal 140 resulting from the 
summation is significantly stronger than the echo difference signal 130 from 
simple pulse inversion. For the case of a linear system, the summation will 
result in a doubling of the scattered power in the echo difference signal. 

In producing a second (inverted) excitation signal 101 5 the phase inversion may 
be applied to the lower frequency conditioning signal, to the higher frequency 
detection signal, or to both. 

An alternative signal processing technique is to obtain a correlation between 
identical first and second signals time shifted relative to one another. 

Figure 15 shows a first dual frequency transmit signal pulse 150 (e.g. of N 
periods of the conditioning signal frequency) and a second dual frequency 
transmit signal pulse 151 which is replica of the first, but shifted by a half 
period of the conditioning signal frequency. 
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Two segments 152, 153 of the same length are then selected from the two 
echoes corresponding to the response to the expansion pressure cycle and to the 
compression pressure cycle. These two segments 152, 153 of the echo signal 
pulses are then con-elated or compared. The resulting correlation is maximal 
5 when no bubbles (non-linear structures) are present (only tissue, or linear 
structures) meaning that the response at the compression and expansion phases 
are similar. The correlation is null (or very small) when bubbles (i.e. non- 
linear structures) are present in the body being examined. 

10 These two implementation options can be performed on the whole received 
echo signals - meaning the raw echo signal is used to perform the summation 
(figure 14) or delay and correlation (figure 15). The signal processing can also 
be carried out on filtered echo signals. A second harmonic filter can be applied 
to the echo signals in either of the figure 14 and figure 15 techniques. 

15 

A super-harmonic filter (third harmonic, fourth harmonic, fifth harmonic or a 
combination thereof) can also be used prior to carrying out the signal 
processing described above. A sub-harmonic filter and or an ultra-harmonic 
filter could also be used. 

20 

Other embodiments are intentionally within the scope of the accompanying 
claims. 



